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Nuclear effects and QCD pcrturbative corrections to the Paschos- Wolfenstein relation- 
*vj ' ship are discussed. We argue that pcrturbative corrections largely cancel out in this 

relationship for total cross sections while the neutron excess correction in heavy nu- 
clei is enhanced by Fermi motion and nuclear binding effects. These observations are 
^ . discussed in the context of NuTeV measurement of the Weinberg mixing angle. 
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^>0 [ The scattering of (anti) neutrino from matter is mediated by charged W^ or W boson 

(charged current, CC), or by neutral Z boson (neutral current, NC). A relation between 
O ■ neutrino-antineutrino asymmetries in the NC and CC deep-inelastic (DIS) cross sections 

r-| ■ was derived long ago by Paschos and Wolfenstein P] 






where 6^ is the Weinberg mixing angle. The derivation of the Paschos-Wolfenstein 
relationship (PW) is solely based on the isospin symmetry and neglects contributions 
from heavy quarks. For this reason this relation is exact for an isoscalar target in a world 
without heavy quarks. In particular, this means that various strong interaction effects, 
including nuclear effects, should cancel out in R~ for an isoscalar target thus making 
Eq. (PQ) a very good tool for the measurement of the mixing angle in neutrino scattering. 

However, in the real world the PW relation is subject to a number of corrections. In 
particular, it must be corrected for the effects due to possible s — s and c — c asymmetries 
in the target (see e.g. Refs. |21E1)- Furthermore, the targets used in neutrino experiments 
are usually heavy nuclei, such as iron in NuTeV experiment j3]. Heavy nuclei typically 
have an excess of neutrons over protons and are non-isoscalar. For a non-isoscalar target 
Eq.(0) is not exact and receives various strong interaction corrections via admixture of 
the isovector components to R~ . In this paper we address QCD perturbative corrections 
and nuclear effects in the PW relationship for non-isoscalar nuclei. 

We will discuss (anti)neutrino DIS in the leading twist (LT) QCD approximation. In 
this approximation the NC and CC structure functions are given in terms of parton distri- 
bution functions (PDFs). In order to simplify discussion of isospin effects, we consider the 
isoscalar, qo{x) = u{x) + d{x), and the isovector, gi(x) = u{x) — d{x), quark distributions 
(for simplicity, we suppress the explicit notation for the Q^ dependence of parton distribu- 
tions). The calculation of the NC and CC cross sections, and the PW ratio in the leading 



order in the strong coupling constant (LO) is straightforward. The next-to-leading order 
(NLO) correction to the PW relation is given in Ref. [2] and the next-to-next-to-leading 
order (NNLO) correction was calculated in Ref. [2]. The result can be written as° 



R~ = --sl, + 6R- 
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where s^ = sin^6'vy, as is the strong coupling, and x~ = J dxx{qa — Qa), with g^ and 
Qa the distribution functions of quarks and antiquarks of type a. The subscripts and 1 
refer to the isoscalar go and isovector gi quark distributions, respectively. In the derivation 
of Eq.Q we expanded in x^ /xq and retained only linear corrections. We also neglected 
contributions due to possible s — s asymmetry discussion of which can be found in Refs. j21 

El El- 

Equation (j21) applies to any, not necessarily isoscalar, nuclear target. We observe that 
SR~ is determined by the valence part of the isovector quark distribution in the target. 
Complex nuclei, such as iron, have unequal number of neutrons (A^) and protons (Z) and 
the isovector quark distribution is finite in such nuclei. In order to understand this effect, 
we first consider a simple approximation which is often used in processing of DIS data. 
In particular, we neglect nuclear effects and view the neutrino scattering off a nucleus 
as incoherent scattering off bound protons and neutrons at rest. We denote qa/x as the 
distribution of quarks of type a in a target T. Then in considered approximation the 
nuclear distribution g^/A is the sum of the corresponding quark distributions for bound 
protons and neutrons 

Qa/A = Zqa/p + Nqa/n- (3) 

We apply Eq.(j2I) to the isovector and isoscalar distributions assuming the isospin invari- 
ance of PDFs in the proton and neutron, i.e. qo/p{x) = qo/n{x) and qi/p{x) = —qi/n{x). 
Then we have qo/A{x) = AqQ/p{x) and qi/A{x) = {Z — N)qi/p{x). For the ratio of average 
light cone momenta in the isovector and isoscalar states, which determine 6R~ in Eq.(P|). 
we obtain 

{xi/xo)a = -5N{xi/xo)p, (4) 

where 6N = {N — Z)/A is fractional excess of neutrons. 

It follows from Eq. (j3)) that the neutron excess correction to R~ is negative for neutron- 
rich targets. In order to estimate the magnitude of this correction for iron target we first 
neglect as terms in Eq.© and compute {x^ /xq)p = 0.45 using the parton distributions 
of Ref.jT] at Q^ = 20 GeV^. Keeping in mind application to NuTeV measurement |3] we 
use 6N = 0.0574 reported by NuTeV |3]. Then we have 5R~ = —0.013. This is a large 
value on the scale of experimental errors of NuTeV measurement since \SR~\ ~ lOcr (for 



"The total cross sections involve the integration of the structure functions over the full phase space 
of X and Q^. Therefore as and the moments x~ of the parton distributions are taken at some average 
scale Q^, which has to be chosen according to specific experimental conditions. The NNLO coefficient 
in Eq.© is given in MS scheme. We also remark that there is an error in the as NLO coefficient in 
Refs. 121 El- I am grateful to K. McFarland for pointing this out. 



a discussion of the neutron excess correction in this context see Ref. jH]). This gives us a 
motivation to study various corrections to R~ in more detail. 

We first discuss perturbative corrections in Eq.Q and consider the difference between 
NNLO and LO approximations 

AR~ = (5i?"(NNL0) - 6R~{L0). (5) 

If we simply use the LO PDFs and as of Ref. [7] at Q^ = 20 GeV^ we obtain from Eq.© 
AR~ = —0.0008, which is about 6% of the LO value of 6R~ . However this calculation is 
not fully consistent, since PDFs as well as the value of a^ in Eg.© should correspond to 
the order of perturbative calculation. It is possible to take into account this effect using 
the results of analysis of Ref. [7j, which provides PDFs to different order up to the NNLO 
approximation. If we do so we observe that the terms in the right side of Eq.(jS)) almost 
cancel each other leading to AR~ ~ 0.7 ■ 10~^. This value is the order of magnitude 
less than the result of a naive calculation (note also that the sign of the correction has 
changed). The reason for this cancellation is that as terms in Eq.Q turned out to be 
balanced by perturbative effects in PDFs which cause the decrease in the ratio xi/xq for 
valence quarks in the proton from 0.457 (LO) to 0.434 (NNLO). 

In order to verify that this cancellation is not accidential we performed similar analysis 
for g2 = 10 and lOOGeVl We have respectively AR' = -1.7-10-^ and 1.9-10"^. These 
values indicate that the cancellation seems to be systematic. It must be also noted that 
such a small value of AR" suggests that the magnitude of perturbative correction to R~ 
is within the variations of 6R~ due to PDF uncertainties of Ref. j7j. Summarizing, we 
conclude that the LO calculation provides a good approximation of R~. 

Now we turn to the discussion of nuclear effects in R~ . In order to improve on 
Eq.Q, we consider nuclear binding and Fermi motion effects (for which we will use the 
abbreviation FMB) in terms of the convolution model of nuclear parton distributions (see, 
e.g., Refs. in 13 CHI)- Then Eq.© should be replaced by 

qa/A = {qa/p)p + {qa/n)^, (6) 

where the two terms in the right side are the quark distributions in bound protons and 
neutrons averaged with the proton and neutron nuclear spectral functions, respectively. 
Similar equation can also be written for antiquark distribution. The explicit expression 
for the averaging in Eq.® is (see piTUj) 

x{qa/p)p = / ded^k Vp{e, ^) ( ^ + ^ ) ^'Qa/pix'), (7) 

x' = ^ = - (R) 

2k -q l + {e + k,)/M' ^^ 

The integration in Eq.Q is taken over the energy and momentum of bound protons (we 
separate the nucleon mass M from the nucleon energy ko = M + e). The quantity Vple, k) 
is the nuclear spectral function which describes the distribution of bound protons over 
the energy and momentum. In Eq.Q, the 2;-axis is chosen in the direction opposite to 
the momentum transfer q = {qo,Oj_,—\q\), and x' is the Bjorken variable of the bound 
proton with four- momentum k. Equation similar to Eq.Q also holds for neutrons with 
the obvious replacement of the spectral function and quark distributions. The spectral 
functions Vp and P„ are normalized to the proton and neutron number, respectively. 



For the isoscalar and isovector nuclear parton distributions we obtain from Eq.® 

go/A = (9o/p)o> (9a) 

?i/A = (^i/p)i' (9b) 

where the averaging is respectively performed with isoscalar and isovector spectral func- 
tions, Vq = Vp + Vn and Vi = Vp — Vn- 

The isoscalar and isovector spectral functions Vq and Vi are very different in complex 
nuclei. In an isoscalar nucleus with equal number of protons and neutrons one generally 
assumes vanishing Vi ^ and nuclear effects are dominated by the isoscalar spectral function. 
In a nuclear mean-field model, in which a nucleus is viewed as Fermi gas of nucleons bound 
to self-consistent mean field, the spectral function can be calculated as 

VMY{e,p)= J2 nx\Mp)\'Sie-ex), (10) 

A<Af 

where (f)\{p) are the wave functions of the single-particle level A in nuclear mean field 
and nx is the number of nucleons on this level. The sum in Eq. (fTUI) runs over occupied 
single-particle levels with energies below the Fermi level Xp- Equation (fTUj) gives a good 
approximation to nuclear spectral function in the vicinity of the Fermi level, where the 
excitation energies of the residual nucleus are small. As separation energy \e\ becomes 
higher, Eq. ljlUp becomes less accurate. High-energy and high-momentum component of 
nuclear spectrum can not be described by the mean-field model and driven by correla- 
tion effects in nuclear ground state as witnessed by numerous studies. We denote this 
contributions to the spectral function as Vcot{£,p)- 

In a generic nucleus the spectral function Vi determines the isovector nucleon distri- 
bution. We now argue that the strength of Vi is peaked about the Fermi surface. It 
is reasonable to assume that Pcor is mainly isoscalar and neglect its contribution to Vi. 
Then Vi is determined by the difference of the proton and neutron mean-field spectral 
functions. If we further neglect small differences between the energy levels of protons and 
neutrons then Vi will be determined by the difference in the level occupation numbers 
nx for protons and neutrons. Because of Pauli principle, an additional particle can join a 
Fermi system only on an unoccupied level. In a complex nucleus all but the Fermi level 
are usually occupied (the Fermi level has a large degeneracy factor). Therefore, Vi is 
determined by the contribution from the Fermi level and we can write 

Vi = iZ-N)\Mp)\'Sie-ep), (11) 

where ep and (pp are the energy and the wave function of the Fermi level. In a nucleus 
with a large number of particles one can use the Fermi gas model to evaluate the wave 
function (pp. In this model \ipp{p)\'^ oc S{pp —p), where pp is the Fermi momentum, and 
we have 

Vi = {Z-N)5{p-pp)S{e-ep)/{47rpl). (12) 



''It must be commented that this statement is violated by a number of effects even in the Z ^ N 
nuclei. The finite difference between the proton and neutron spectral functions is generated by Coulomb 
interaction and isospin-dependent effects in the nucleon-nucleon interaction. The discussion of these 
effects goes beyond the scope of this paper and we leave this topic for future studies. 



We now apply these equations to calculate the binding and momentum distribution 
effects on average quark light-cone momenta in the isoscalar and isovector quark distri- 
butions. Integrating Eq.((7j) over x and keeping the terms to order e/M and h^/M"^ we 
have 

^ = -.7v(l + rf£).v,. (13b) 

where £0 and Tq are the separation and kinetic energy per nucleon averaged with the 
isoscalar nuclear spectral function Vq and Tp = p|./(2M). In order to quantitatively 
estimate this effect we observe that the energy of the Fermi level ep equals the minimum 
nucleon separation energy. For the iron nucleus we take ep = — lOMeV and pp = 260 MeV 
(the corresponding energy Tp = 36 MeV). In order to calculate the isoscalar parameters 
Eq and To we use the model spectral function of Ref . ^J which takes into account both the 
mean-field and correlated contributions (see also 521)- We find that the naive neutron 
excess correction by Eq.(j3]) should be increased by the factor 1.055. 

We now discuss these results in the context of NuTeV effect jl]. We assume that 
the weak mixing angle can be calculated from Eq.Q in terms of experimental R~. In 
particular, we are interested in the variation of s^ because of nuclear effects and effects 
of higher order in ag, since NuTeV analysis was carried out in LO approximation without 
nuclear effects. The correction As^ is apparently given by the difference between the 
corrected and uncorrected experessions for 6R~ 

As^ = (5i?-(NNL0+FMB) - 6R-(L0). (14) 

Perturbative corrections are largely canceled out in the difference as discussed above and 
the resulting value As^ = —0.00065 is practically saturated by nuclear effects. With 
certain care this value can be viewed as a correction to the value of the Weinberg angle 
sin^ 9\Y measured by NuTeV jlj ? 

In summary, we discussed perturbative QCD corrections to the PW relationship to- 
gether with nuclear binding and Fermi motion effects. A cancellation of QCD perturba- 
tive corrections to the PW relationship for the total cross sections has been observed. We 
found a negative correction due to nuclear effects to the PW relationship for the total 
cross sections for neutron-rich targets and estimated this effect on the Weinberg angle of 
NuTeV measurement. 
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■^It should be remarked in this context that total cross sections were not measured in NuTeV experi- 
ment. For this reason a more detailed analysis of differential cross sections with the proper treatment of 
experimental cuts and (anti)neutrino flux is needed. 
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